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THE GAS PHASE IN POROUS FUEL-CELL ELECTRODES 


by Robert W. Easter 
Lewis Research Center 

SUMMARY 


A generalized physical model of a porous fuel-cell electrode has been presented 
which does not require postulation of any particular geometrical configuration. This 
allows the writing of the basic conservation equations describing the electrode in rel- 
atively simple general forms. 

Applying the equations in a number of situations shows that the composition of the 
gas phase in fuel- cell electrodes is dependent on the values of the three parameters J,G, 
and (p, where 


J = 


"Ivol^ 


nRfCTDn^O-R^ 6A 


D 


6 =■ 


I-R 


D 


; and 


I-HgO 


D 


cp =. 


I-R 


^HgO-R 

T is electrode tortuosity, I^^j^ is total current produced by electrode volume element, L 
is z- value of greatest penetration of gas phase toward bulk electrolyte, is equiva- 
lents per gram-mole or reactant, f is Farraday's constant (9.65x10“^ A-sec/g-equiv), 
C^ is total concentration, e is ratio of cross-sectional area to geometric cross- 
sectional area, and 6A is increment of electrode surface area. The quantity J may be 
thought of as a dimensionless current density, while G and 0 are ratios of binary-gas- 
phase diffusion coefficients. The significance of these parameters is that for many pur- 
poses they make detailed analysis unnecessary. Considerable information is given sim- 
ply by the m^nitude of the parameters. 

In particular, criteria are established whereby it can be determined whether gas- 
phase concentration gradients within the electrodes are significant. The parameters are 
also related to the accumulation of gaseous inert impurities in the electrode. 



INTRODUCTION 


Vital to the operation of porous fuel-cell electrodes are large amounts of gas-liquid 
and liquid-solid interfacial area in close proximity to each other. Thus, understanding 
of fuel-cell operation requires knowledge of the processes going on in all three phases: 
gas, liquid, and solid. Because of the interaction of these processes, any description of 
fuel- cell operation is of necessity quite complicated. A further hindrance to under- 
standing is the fact that only overall performance may be investigated experimentally for 
most electrodes. 

Although a number of fairly complex mathematical models of the operation of porous 
fuel- cell electrodes exist (see refs. 1 and 2 for summaries) none of these consider the 
gas phase within the electrode in any detail. This neglect is due chiefly to the complex- 
ity of porous electrode behavior. It is difficult enough to describe kinetic and liquid- 
phase phenomena, without coupling possible gas-phase effects. 

In this work, however, an attempt is made to shed some light on the state ol the gas 
phase in operating electrodes. In order to focus on the gas phase, and to avoid com- 
plexity, the following approach was taken: 

(1) Establish a generalized physical picture which can be analysed in an orderly 
manner. 

(2) Present the basic gas-phase conservation equations which apply to the established 
picture. 

(3) Obtain the specific forms of these equations that apply to specific electrode types. 

(4) Solve the specific equations and draw suitable general conclusions from the re- 
sults. 

Rather than solve simultaneously the analogous equations describing the liquid and 
solid phases in the electrode, reasonable forms for the liquid- and solid-phase phenom- 
ena are postulated. 

While an approach of this nature usually limits the confidence with which results may 
be applied in cases which differ from those considered, it will be shown that several 
characteristics appear consistently, regardless of the specifics of a particular case. 
Furthermore, the first two sections of this report, corresponding to points (1) and (2) 
in the preceding list, are quite general and form a useful basis for analyzing the gas 
phase in any electrode. 

There are several ways in which the information generated by analysis of the gas 
phase is useful: 

(1) In the construction of complete models of electrode behavior, simplifying assump- 
tions regarding the gas phase may be made more intelligently. It is possible to deter- 
mine, for instance, the conditions under which concentration gradients are approximately 
linear, or negligible. 
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(2) The consequences of impure fuels and oxidant may be assessed. 

(3) The effects of the bulk gas in the gas cavity on the electrode interior are more 
recognizable. 


GENERALIZED PHYSICAL MODEL OF ELECTRODE INTERIOR 

To analyse the mass transfer characteristics of a porous fuel- cell electrode it is 
necessary to establish an idealized picture of the electrode interior. If the results of 
the analysis are to have wide application, the picture or model adopted should be as gen- 
eral as possible. Thus, the parameters to which the model gives rise, while being few 
in number, should have meaning regardless of the specific electrode under considera- 
tion. To this end, consider a cylindrical section "punched out" of an operating elec- 
trode in a fuel cell, with a diameter large with respect to the mean diameter of the par- 
ticles of which it is constructed, but small with respect to the overall electrode meas- 
urements (see fig. 1). The volume of this cylinder will be filled with solid material 
(catalyst; support structure), liquid, and gas. Current is produced at sites on the solid 
surfaces which (1) are catal3dic in nature, (2) are immersed in liquid which is continu- 
ously connected by a liquid path to liquid covering a similar site in the other electrode, 
(3) are connected by an electronic path to active sites in the other electrode, (4) are rel- 
atively near an area of gas-liquid interface. The liquid will consist of a solvent which 
may take part in the electrode reaction, charged species which may appear or disappear 
at reaction sites and dissolved gases. The gas phase will be made up of reactants sup- 
plied externally, vapor of the volatile liquid species, and possibly other species present 
as impurities in the reactant supplied. 

The bulk or net movement of the reactant gas in the electrode will be in a direction 
away from the gas cavity and normal to the electrode face. The twisting and turning 
nature of that movement (and of the movement of other gas -phase species) may be ac- 
counted for by specifying the "tortuosity" of the electrode, whereby the effective path 
length for transport is the product of the geometric distance and the tortuosity. The 
(perhaps position dependent) ratio of gas-filled cross-sectional area to geometric cross- 
sectional area e(z) defines the relative "openness" of the electrode. 

Transfer of mass between phases takes place at areas of gas-liquid interface. If 
current is being produced in the vicinity of a portion of interfacial area, there will be a 
net transfer of reactant gas across the interface at that area. If no current is being pro- 
duced nearby, reactant equilibrium will exist between the two phases at that area, in the 
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sense that there will be no net transfer of reactant across the interface. In either case 
there may be net transfer of some other species across the interface. 

The foregoing discussion leads logically to the definition of a current density based 
on a unit area of gas-liquid interface. The distribution of reactant gas in the gas phase 
can then be determined if this current per imit interfacial area can be specified and if 
the’ occurrence of gas and liquid phases in the volume under scrutiny can be described. 

Expressions for the interfacial current density may be derived, but require some 
additional specification of electrode structure, reaction kinetics, and liquid behavior and 
so lose some generality. Considerable information can be had, however, without fixing 
the detailed interfacial current density distribution, but rather by investigating various 
postulatory forms of the distribution. Generally, then, this work will assume various 
forms for i^(z), where i^ is the interfacial current density and z the coordinate dis- 
tance from the electrode surface. 

Finally, we shall define an important quantity, ^(z), where ^(z) is the position- 
dependent quantity: interfacial area per unit geometric volume. It is this parameter 
that preserves the generality of this approach. It has significance regardless of whether 
it is considered an unspecified parameter or whether it is applied to a specific idealized 
picture of an electrode (e.g. , the simple pore model (ref. 3), or the flooded agglomer- 
ate model (ref. 4)). 

The situation within a "punched out" volume of electrode, with particular regard to 
the gas phase, can now be analyzed. Note that the molar reaction rate per unit volume 
is [ij^(z)/npjf] |(z) where Up^ is equivalents per mole, and that the total current produced 
by the volume is 

^vol" I ii(z)?(z) 6Adz 

where L is the depth of electrode in which there exists some gas phase and 6 A is tne 
geometric cross-sectional area of the cylindrical volume. (All symbols are defined in 
appendix A.) 


EQUATIONS DESCRIBING BEHAVIOR 

To completely describe mathematically the transport processes in the gas phase one 
needs (1) a continuity or conservation equation for each species, (2) expressions relating 
flux and concentration of each species, (3) an equation of state for the mixture, (4) an 
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equation of motion or momentum conservation, and (5) an equation of conservation of 
energy. 

The general form of the continuity equations to be used herein is 


-a(Nie) 

r?.(z)|(z) 

0z ^ 


e SC. 


St 


( 1 ) 


where is the z-direction flux, based on gas-filled area, of species i with respect 
to fixed coordinates; is the rate in moles per unit area of interface at which species 
i is transferred into the liquid phase; and is the gas-phase concentration of species 
i. This expression embodies the assumption that the net transfer of material out of the 
volume through its sides is negligible; that is, reactants and products are exchanged 
with the surroundings only in the positive or negative z-directions. 

Equation (1) is then an effectively one-dimensional expression; variations from re- 
gion to region in the electrode may be accounted for by suitable boundary conditions at 
the electrode surfaces. 

The fluxes and concentrations of the various species are related by the so-called 
Stefan-Maxwell equations (ref. 5, p. 570). These consist of an equation of the following 
form written for each of the n species: 



T SZ 



^T^ij 


(X.N. - X.N.) 


( 2 ) 


where Xj is the mole fraction of the component, C,p is the total concentration, t 
is the electrode tortuosity described earlier, and the diffusivity of the pair i-j in 
a binary mixture. This expression embodies the assumptions that the gas phase is an 
ideal-gas mixture and that transport takes place by ordinary diffusion (and possibly con- 
vection) only. The utility of the Stefan-Maxwell equations lies in the fact that they ex- 
press the fluxes in multicomponent systems in terms of a quantity (C-pD-j) that may be 
calculated from binary system data or correlations. Furthermore this quantity is inde- 
pendent of both pressure and relative compositicn for ideal gases at low pressures (ref. 

5, p. 510). Note that the equations reduce to the usual form for the binary situation, and 
that the use of the Stefan-Maxwell equations requires that the equation of state be the 
ideal- gas law. 

The equation of motion may be written in the following form, which is known as 
Darcy's Law. Gravitational effects are neglected. 
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( 3 ) 


k ^ 
II dz 


where is the so-called superficial velocity - the volume rate of flow through a unit 
cross-sectional area; p. is the viscosity of the mixture; and k is the permeability of 
the porous material. The volume rate of flow per unit gas -filled area is given by 1/e 
times equation (3). Equations (1) and (2) may be solved for the gas-phase compositions 
without need of equation (3) if the species fluxes are known. This is the approach to 
be taken in the remainder of this work and thus equation (3) will not be dealt with further. 
In a rigorous analysis, however, equation (3) would be solved simultaneously with equa- 
tions (1) and (2) to give the net convective mass flow. 

An equation of energy in steady-state form may be written for the gas -filled portion 
of the control volume: 


dz 



+ ?(z)Q(z) = 0 (4) 


where is the thermal conductivity of the gas mixture, T^ the local gas temperature, 
hj^ the molar enthalpy of the species i, and Q(z) the rate at which energy is transferred 
across a unit area of gas-liquid interface. This expression assumes that kinetic and 
potential energy effects are negligible and that the z-direction energy flux is made up of 
ordinary conductive and convective components only. 

This work will not be concerned with solving expression (4). As a rule, tempera- 
ture gradients do not directly affect mass transport, althrough the converse is not true, 
as evidenced by the appearance of the species fluxes in equation (4). Equations (1) 
to (3) are written to show temperature dependence only in the physical properties, and 
these can be evaluated at the mean or nominal temperatures chosen without too much 
error. 

A further assumption is made regarding these descriptive expressions in the re- 
mainder of this work - all time-dependent terms are negligible. This means that the 
example solutions presented cannot describe situations that are clearly transient (e.g. , 
behavior immediately after a change in current drawn). They can, however, be applied 
to "slow” processes, such as the gradual buildup of inert species in the gas phase, by 
considering the process to proceed as a succession of pseudo-steady-states. 

Finally, in all the specific cases considered in the following sections the quantity 
^(z) is assumed to be constant over the range z = 0 to z = L. In physical terms this 
means that the cross-sectional area available for gas-phase transport is constant 
throughout the region of the electrode in which the gas phase is present. 
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If the liquid in this region is confined to small pores within the solid particles 
and/or thin films upon the particles, then e may be identified with the so-called 
macroporosity of the electrode. If, however, considerable amounts of liquid are pre- 
sent in partially liquid-filled large interparticle pores and consequent menisci, the 
available area for gas transfer will be less than that given by the electrode macropor- 
osity, since that property is based on dry electrode structure. 

Variations in liquid distribution may be investigated if desired by postulating forms 
for the z-dependency of e and solving the pertinent equations by numerical methods. 

The remainder of this work is concerned with the analysis of various operating sit- 
uations using the expressions developed in this section. The major goals in each case 
are (1) to assess the validity of the assumption of constant reactant partial pressure 
throughout the gas phase, and (2) to investigate the distribution of inerts (if present) in 
the gas phase in the electrode. 


ANALYSIS OF SPECIFIC SITUATIONS - THE ANODE AT WHICH 
WATER IS BEING REMOVED 

Here the gas phase is assumed to be made up of hydrogen gas, water vapor, and an 
inert impurity. The composition of the gas in the gas cavity immediately above the con- 
trol volume is assumed to be known. It should be borne in mind that the gas cavity com- 
position may differ from place to place and that the analysis of a single volume incre- 
ment may not represent the conditions over an entire electrode. 

For hydrogen, integration of equation (1) yields 


N 


H 


^vol 
2 2f6A€ 


/ z 

Vtj ^ dz 
2 


Recall that the molar reaction rate of hydrogen perunit volume is equal to i^^/2f. 

At steady state this reaction rate per unit volume may be set equal to the rate of inter- 
phase transfer of H2 per unit volume. The two rates will be unequal only if there ex- 
ists considerable z-direction transport of dissolved hydrogen in the liquid phase, an un- 
likely circumstance because of the very low solubility-diffusivity product of Hg in elec- 
trolyte solutions. 

Thus rju = i;/2f . And from the preceding equation 
^2 ^ 
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(5) 



^vol 
2f 6A e 



dz 


If water vapor is being transported in the gas phase within the electrode, equation 
(ly may be integrated to yield 




vol 
6A e 


(6) 


If the impurity solubility is assumed to be negligible, the impurity continuity equa- 
tion (eq. (1)) becomes simply Nj = 0. 

Taking advantage of the fact that Xjj q = 1 - Xjj - Xj, and since Nj = 0, the Stefan - 

2 2 

Maxwell equations may be written (from eq. (2)): 


1 dXj 


Xj dz 


N„ +0N„ 
^2 ^2 




(7) 


and 


dX, 


dz 


-- X 


H 




%2 * %0 


^T^I-H. 


N [.y.Xj(l-v)] 


( 8 ) 


where 


D 


I-Hn 


e =. 


D 


I-HgO 


(p = 

%0-H2 

Now, if the distributions of current density, interphase water transfer, relative gas- 
filled cross-sectional area, and active interfacial area are known (along with the re- 
quired physical properties), the integrals in equations (5) and (6) may be evaluated. 
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equations (5) and (6) substituted in equations (7) and (8), and equations (7) and (8) solved 

for Xj(z) and (z). 

2 


K the distributions are unknown, valuable information can be gained by postulating 


various forms of i. , r?, 


HoO’ 


and ^ which appear capable of describing the situation in 


a particular electrode or under a particular operating condition. In the following pages 
three cases of this sort are considered. The relation between each case and actual 
operating modes will be discussed and the general conclusions which may be drawn 
will be pointed out. 


Case A 

Essentially all the current is produced at or very near the bottom of the control 
volume. Likewise, all interphase transfer of H2O takes place approximately at the 
bottom of the volume. 

The gas-phase transport in the volume is then an equimolar countercurrent diffu- 
sion process. The situation described here is approached in nonwetting electrodes or in 
partially wetting electrodes operating in a "dried-out" condition. 

Mathematically, 


N. 


H 


Vol 
2 2f 5A e 


= - N. 


HgO 


and thus for a "pure" gas phase (i. e. , Xj = 0), equation (8) yields 


Xh_(z) 




( 9 ) 


where X^^ is the mole fraction of hydrogen in the gas cavity (at z = 0) and 


J = 


"'vol^ 


2fC,pDjj^Q_jj^e 6A 


( 10 ) 


This quantity J will appear in all cases (as will 6 and (p if impurities are pres- 
ent) and is of prime importance in characterizing the gas phase. Note that the quantity 
I^oi/SA is equivalent to the apparent electrode current density. 
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In this case the composition change with z is directly proportional to the magni- 
tude of J. For example, the composition change from top to bottom of the volume is 
(from eq. (9)) 



( 11 ) 


If the gas phase contains inert impurities, their distribution may be had by solving 
equation (7), yielding 


Xj(z) = X|^ exp 




(12) 


The substitution of this expression into equation (8) has the effect of adding another 
term to the expression for the reactant mole fraction, that is, now 


Also, 

% - = J * - '} 

Values of 9 and <p at two temperatures are given in table I for four likely impuri- 
ties. For ideal-gas mixtures at low pressures, 9 and (p are independent of composi- 
tion and total pressure (see appendix B for calculation of diffusivities). 

In case A, it was postulated that current production takes place at the bottom of the 
volume. In this situation the mole fraction (or partial pressure) of H 2 at z=L should 
be as high as possible. The higher the partial pressure the greater the concentration of 
dissolved hydrogen in the liquid phase, and thus the lower the electrode polarization for 
a given current density. 

The partial pressure of water vapor at the bottom of the volume is also of interest, 
since it is related to the concentration of the electrolyte "beneath" it. For example, if 
the gas cavity composition is fixed, at steady-state the partial pressure of H 2 O at z=L 
must be that which allows vapor-phase transport of the stoichiometric amount of H 2 O 
out of the electrode. And this dictates what the concentration of the liquid electrolyte 
immediately below must be (for a given temperature) . 

Evaluation of equation (12) at z=L and equation (14), together with the relation 
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Xj +Xjj Q + Xjj = 1, yields the composition of the gas phase immediately above the 
2 2 

electrolyte. Figures 2(a) to (c) give example results for the four impurities in table I. 
The temperature was chosen to be 343 K and the HgO vapor mole fraction in the gas 
cavity above the volume is constant at 0. 15. The J = 0 points correspond to open- 
circuit, static equilibrium situations and thus represent the constant gas cavity com- 
position associated with each curve. Note that figure 2(b) holds for both N 2 and Ar 
and that all three figures are for J ranging from 0 to 0.2. 

Although the three figures are rather different in appearance, one important simil- 
arity exists - for a given gas cavity composition and J the value of X„ is nearly the 

"2 

same regardless of the impurity species. This is because of the relative ease with 
which Hg diffuses through gas mixtures. The binary diffusion coefficients of Hg with 
another gaseous species are relatively high and relatively insensitive to the identity of 
the other gas. This effect is shown clearly in figure 3 wherein X„ (z=L) is plotted as 

"2 

a function of J for the several values of X^^. Here Xjj (L) is essentially identical 

2 


for CH^, Ng, and Ar, the exception being He, which shares with H 2 the ability to dif- 
fuse with relative ease. 

Although the mole fraction of H 2 in the gas phase at z=L is more or less inde- 
pendent of impurity species, the mole fraction of H 2 O is not. In fact, the highly 
species-dependent balance between H 2 O and impurity accounts for the marked differ- 
ence in appearance of the three parts of figure 2. In general, the higher the value of 9 
(see table I) the greater the concentration of water vapor at z=L for a given gas cavity 
composition and J-value. This is an important relation because it gives information 
about the relative dilution or concentration required of the electrolyte at a given steady- 
state condition. Thus, if methane is introduced into a formerly pure gas phase where 
the partial pressure of HgO in the gas cavity is controlled, the electrolyte must be con- 
siderably diluted if the same amount of current is to be produced (and product water re- 
moved) on a steady-state basis. But introduction of helium, on the other hand, actually 
leads to a slightly more concentrated liquid electrolyte. 

The results of this case have been discussed in some detail - not because of its 
particular importance but rather to illustrate the sort of information yielded by solutions 
of the describing equations. The following cases are considered in somewhat less detail. 


Case B 

The anode of case A was operating quite inefficiently from the standpoint of catalyst 
utilization. Now consider the other extreme - a cell where current is produced through- 
out the control volume at a constant rate per unit volume. The transfer of H 2 O into the 


11 



1 II I 


gas phase will be assumed to take place at the same rate and this can be seen to be a 
situation where transport in the gas phase is again through equimolar countercurrent 
diffusion. A cell behaving in this manner would have to have (1) an electrolyte with very 
high conductivity, (2) very good liquid- phase mass transport characteristics, and (3) 
large cross-sectional area liquid paths between points in the anode and points in the bulk 
liquid. 

Mathematically, for this case 


-i- I = - rjjj = Constant 
2f 2 


^vol 
2f 6A L 


and so from equations (5) and (6) 



(i 

2f 6A e V 



For a "pure" gas phase, solution of equation (8) yields 


X„ (z) 
^2 





and now 


- X. 


"H, 


He 



(15) 


The distribution of any impurities present is now (from solution of eq. (7)) 



(16) 


and substitution of this into equation (8) yields the solution for H2 distribution in the 
presence of inerts: 




J 1 - 




(17) 
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Equations (15) to (17) can be seen to be identical to equations (9), (12), and (13), ex- 
cept that where J(z/L) appeared previously, the quantity 

L 2 \L/ 


is now found. The analog of equation (14) is now 



Now since current is being produced throughout the volume, the composition at z=L 
is no longer of prime importance. Comparison of the corresponding figures, however, 
shows that the composition change from z=0 to z=L will be only about half as great in 
this case, compared to case A, for the same gas cavity composition and J value. Fig- 
ure 4 shows the composition at z=L from equations (16) and (18) for Ng. Compare 
this figure to figure 2(b). 

In this case the integrated mean composition is more indicative of the electrochem- 
ical behavior of the electrode than is the composition at z=L. From equation (16) 






dz 


This expression may be rewritten 


X 


I 





(19) 


where u = z/L and xjy = - 0)J/2. 

The expression (from eq. (17)) for the mean mole fraction of H 2 becomes 




- J + X 
3 


gc 

I 




and the mean mole fraction of H 2 O 


is given by 


( 20 ) 
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( 21 ) 


1 - Xj - Xh^ 

Equation (19) can be easily integrated numerically to give Xj, and the mean compo- 
sition of the gas phase can then be found. Figures 5(a) to (c) give the mean composi- 
tions of the gas phase under operating conditions paralleling case A. 

By comparing figures 2 and 5, it can be seen in general that for a given impurity 
the composition in case B is equivalent to the composition at a lower J-value in case A. 
In other words, the curves for a given impurity in case B are essentially identical with 
the lower portions of the case A curves. Therefore, the remarks made earlier regard- 
ing the dilution and concentration of the electrolyte apply in this case, although now it is 
the mean concentration of all the liquid in the electrode about which the figures yield 
information. 

Before extracting any generally applicable conclusions it would be well to investi- 
gate a case where the gas -phase diffusion is not equimolar and countercurrent through- 
out the electrode. Deviations from equimolar countercurrent behavior occur if the rate 
of evaporation of H 2 O per unit volume is different from the rate of current production 
per unit volume. 


Case C 


Let the current be produced at a constant rate throughout the control volume, as in 
case B. Now, however, it will be assumed that all the interphase transport of H 2 O 
takes place at or very near the z = 0 plane. Thus, the hydrogen now diffuses through 
a quiescent mixture of H 2 O (and inert). This is an extreme case but it serves to bracket 
real situations and this condition might be approached in a cell under certain circum- 
stances. This anode would necessarily have the properties described in case B and 
would, in addition, require that the H 2 O produced in the upper half of the volume dif- 
fuse through the electrolyte to the z = 0 surface before evaporating. (This is the sim- 
plest way of meeting the conditions of this case while satisfying the steady- state stoichio- 
metric requirement that one-half the water produced at the anode be transported via the 
liquid phase to the cathode while the other half is removed continuously. ) 

Mathematically, ’’In ^ Njj q ~ that equation (1) pertains to this 

2 2 

case rather than equation (6) since there is no H 2 O vapor transport within the electrode. 

As in case B 


^ = Constant = — — 
2f 2f5A L 
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and 



(i 

2f 6Ae V 



This yields from equation (7) 


Xj - Xp exp ^ 


---(-V 

L 2 VL/ 


( 22 ) 


and the expression for the H 2 distribution from equation (8) is 


X 


He 


1 - xf ^ exp <1 ^ 

^ ’ cp 


L 2l,L/J 




L 2\L/ 


(23) 


This expression also is applicable when the gas phase is pure, that is, when X^^ = 0. 

Figure 6 shows the composition at z = L with N 2 as the impurity. This figure 
may be compared with figures 2(b) and 4. In this case the mole percent of water vapor 
is very nearly constant, varying only very slightly with position and J- value. 

Equations (22) and (23) may be integrated numerically to get the mean composition of 
the gas phase. Figure 7 shows the result for the, by now familiar, example conditions. 
The behavior is approximately the same for CH^, N 2 , or Ar since their respective 
values are nearly the same. Again, the mole fraction H 2 O is nearly unchanging. 

The mean composition when He is the impurity varies only a few mole percent from 
J = 0 to J = 0.2 and would be superimposed on the low-J-value portions of the curves 
of figure 7. 

In contrast to cases A and B, the mole fraction of impurity now increases as J in- 
creases. This is to be expected, however, because now the net flux of mass in the pore 
is proportional to the current density and thus increasing amounts of inerts are carried 
into the pore with increased reactant fluxes. 


Discussion of Anode Results 

Three modes of operation have been investigated for a porous fuel- cell anode with 
hydrogen fuel and an aqueous electrolyte. These three cases represent extremes which 
bracket the situations found in most existing anodes. 


15 



Most practical anodes are completely or partially hydrophilic and thus liquid may 
be distributed fairly evenly throughout the electrode. Except possibly under optimum 
conditions, however, there are appreciable electrical potential gradients in the 
z-direction. And thus while current generation is not confined to the very bottom of the 
electrode, neither is it likely to be distributed evenly throughout the electrode. 

Although the rate at which H 2 O is evaporated across an area of gas -liquid inter- 
face should be proportional to the amount of current generated in the vicinity, the ratio 
of interphase H 2 O transfer to current density is probably somewhat greater at low 
z-values than at values of z near L. 

It should also be recognized that in an actual fuel cell the pattern of gas flow in the 
gas cavity above the anode may be such as to force different portions of the anode into 
different modes of operation. If, for instance, a dry stream of reactant is used to re- 
move product water vapor, the portion of electrode near the gas cavity inlet may be con- 
siderably drier than the portion near the exit. The region near the inlet may thus be 
more nearly represented by case A, while the region near the exit has perhaps more the 
character of case B. Also, if the gas flow is clearly directional for some length the gas 
cavity composition will vary along this length, becoming depleted of reactant in the di- 
rection of flow. 

It is clear then that a rigorous mathematical description of a fuel- cell anode in op- 
eration requires simultaneous solution of kinetics, gas- and liquid- phase transport, and 
gas cavity flow. The idealized cases A to C offer an alternate means of gaining gener- 
ally applicable conclusions, however, by virtue of certain similarities between them. 

In each case, expressions have been derived for a gas-phase composition which is 
representative of the electrode. In case A it is the composition immediately above the 
gas-liquid interface, and in cases B and C the integrated mean composition in the elec- 
trode. It can be seen that in each case the same three parameters J, 6 , and <p appear 
and that their effects are the same. 

Thus, regardless of the detailed mode of operation of an anode, the state of the gas 
phase within the anode and the effect of this gas phase on performance can be character- 
ized if the values of J, 0, and ip are known. 

The following conclusions apply. They are first listed, and then discussed at length 
below ; 

(1) The larger the value of J the more position-dependent is the partial pressure of 
reactant. 

(2) If inert impurities are present, the higher the value of 0, the greater the 
position- dependency of the water vapor partial pressure. 

(3) If inerts are present, the higher the value of <p the more likely is the impurity 
to be distributed equally throughout the gas phase. This conclusion is related to con- 
clusion (2) and thus will not be discussed separately. 
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With regard to conclusion (1), it can be seen from the figures that case A shows the 
greatest variations in composition and so it may be used as a limiting case in testing 
values of J. 

More specifically, if J is less than 0.05, the mole percent of H 2 should be within 
6 mole percent of the gas cavity value everywhere in the electrode regardless of inerts. 
The range 0 s j < 0. 05 does in fact encompass most practical electrodes operating at 
current densities of less than 1 ampere per square centimeter, and so in most cases the 
assumption of constant reactant partial pressure in the anode is a good one. High cur- 
rent density, relatively closed structures (high value of r/e), or low temperatures 
(which decrease Dtt however, increase J. 


Conclusions (2) and (3) also deserve a more detailed discussion. If the gas phase is 
pure, it is obvious that the decrease in mole percent H 2 from the gas cavity value will 
equal the increase in mole percent of HgO vapor. The situation is apparently more 
complicated if inerts are present. Although in every case the characteristic composi- 
tion shows an increase in water vapor with increasing J, the magnitude of the increase 
(from the gas cavity composition) is highly dependent on 0, that is, the impurity species. 
With the exception of case C, the concentration gradient of water vapor at steady-state is 
determined by the amount of water vapor which must be transported from the electrode 
interior to the gas cavity and by the relative ease with which this transport can take 
place. A large value of 9, then, may be interpreted as a propensity to impede HgO 
vapor diffusion, forcing the formation of a relatively larger partial pressure gradient to 
transport the stoichiometrically correct amount of water vapor. 

It is difficult to provide quantitative guidelines for 0-values since even CH^ (the 
largest 0 of the species in the examples) has little effect on a cell operating in the mode 
of case C. Thus, unless an anode is known to operate in a particular manner, the best 
course is to assume its behavior is approximated by case A, and find the representative 
composition from figures 2 and 3. Since case A represents a worst case it yields an 
upper limit on the dilution effect regardless of the anode operating mode. 


ANALYSIS OF SPECIFIC SITUATIONS - THE CATHODE 
WITHOUT H 2 O REMOVAL 

Here the gas phase is assumed to consist of O 2 gas, H 2 O vapor, and an inert im- 
purity, Again, the gas cavity composition is assumed to be known. In some situations, 
however, it might be preferred to specify the H 2 O partial pressure at z=L, (i.e., the 
vapor pressure of the bulk electrolyte), and the mole fraction of inerts in the gas cavity. 
The equations may be solved with these boundary conditions as easily as for a known gas 
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cavity composition, and herein the gas cavity composition is specified merely for con- 
venience. 

The continuity equations (eq. (1)) take the forms 



-JLf i I dz 
4f 6A e MejQ ^ 


N 




dz 


(24) 


(25) 


Note that 77 


H2O 


is a positive quantity if 


H2O is entering the liquid phase. 


Again using the pseudo-steady-state assumption, Nj = 0. 

The chief difference between these expressions and the corresponding expressions 
for the anode is that now there is no net removal of H2O from the volume via the gas 
phase. Thus, equation (25) is concerned with vapor-phase transport of H2O from one 
part of the volume to another. 

The Stefan-Maxwell equations (eq. (2)) take forms analogous to the corresponding 
forms for the anode (compare with eqs. (7) and (8)): 


i_^^ = 1 /Nq 

Xj dz C^Dj q V ^2 "2^ 


(26) 


and 


dX, 


- X, 


dz 2 


No,) = - TW— ^0, + X,(l - *.)] (27) 


where now 



^1-H20 


and 
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cp = 



%2O-02 


Here again the procedure is to postulate various spatial distributions of interphase 
transfer and to investigate the results with the aim of drawing generally applicable con- 
clusions. 


Case A 

Essentially all the current is produced at or very near the bottom of the volume. 

The reactant gas thus is transported through a stagnant mixture of H2O vapor and inert. 
This case represents a hydrophobic, or a dried-out, electrode. Thus equations (24) and 
(25) become 



Vol 

4f 6A e 


and 


NH 2 O = 0 


Substituting these expressions into equation (26) yields the inert distribution: 




where 


J = 


^IvolL 


4CtDh 2O-02^ 6A e 


The J here defined is similar to the J defined for the anode with 4C,pDjj q_q 

2 2 

now appearing instead of 2CrpDjj q_jj . Because Djj q_q is only about 1/4 Djj q_jj 

2 2 2 2 2 2 
for a given temperature, J for the cathode is about twice the magnitude of J for the 

anode for identical geometric factors and current densities. 
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Equation (27) may be solved for the reactant distribution: 

Xo^(2) = 1 - Xf <= exp |£ - 4=0 exp |j (28) 

Note that since there is no transport of water, vapor, the parameter 9 does not 
appear. 

The most interesting feature here is that both the water vapor and the inert tend to 
be present at higher concentrations at z = L than at z = 0. For the inert in particular. 


Xj(L) 

Xf^ 



(29) 


Values of cp and 6 for the same impurities and temperatures as the anode values 
are given in table 11. Figure 8 shows Xj(L)/X^^ for the impurities in table n at 343 K. 
This "pile up" of inerts at the bottom of the volume results from the net mass flow (of 
reactant) in the z-direction, which carries along some of the inert. 


Case B 


The next case for consideration, as a contrast to case A, is the situation where the 
current is produced at a constant rate per unit volume for all z. As in the preceding 
case there is no HgO vapor transport. The mathematical expressions for this case are 
identical with those in case C for the anode. Thus, equations (22) and (23) describe the 
gas-phase composition in this case for O 2 as the reactant as well as for H 2 . Equation 
(22) evaluated at z = L now yields 


Xj(L) 




(30) 


And so, although the inert still tends to accumulate in the bottom of the volume, the 
ratio of equation (30) is equal only to the square root of the ratio in equation (29). 
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Case C 


Here a somewhat closer approximation to reality is to let the current production in- 
crease linearly with z. Thus 




2 ^ 
6A L L 


Again, assuming no H 2 O vapor transport, N 


HgO 


= 0 


and from equation (24) 



^vol 
4f 6Ae 



The tendency for the majority of the current production to be near the bottom of the 
volume should be more pronounced in the cathode than in the anode because of reactant 
solubility gradients. In the anode the electrolyte in the upper region (z < L/2) is dilute 
relative to the lower region and so the gas solubility is relatively higher in the upper re- 
gions. This effect, which would tend to produce more evenly distributed current produc- 
tion, is reversed in the cathode; so the potential gradient, the availability of liquid spe- 
cies, and reactant solubility all tend to favor the lower region of the cathode. 

The inert distribution from equation (26) is now 


Xj - Xp exp 




and the reactant distribution (from eq. (27)) 




z 

L 







(31) 


(32) 


As before it is desirable to specify the gas-phase composition at some point which 
is representative of the electrode as a whole. In this case, since current production is 
greater near the bulk electrolyte than near the bulk gas (i.e. , greater for larger z/L), 
this composition should be a weighted rather than a simple mean. To this end the com- 
position at z/L = \/ 1/2 was chosen, which represents the point at which one-half the 
total current has been produced. Figure 9 shows this composition for Ng as the inert 
at 343 K. The mole fraction of HgO vapor in the gas cavity is held at 0. 10. Note that. 
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because of the similarity of <;o-values, this also represents (to within a few percent) the 
situation when CH^ or Ar is the impurity. 

It is also instructive to use the results of this case to compare a ’’pure" oxygen 
cathode to an air cathode. Figure 10 shows the mole fraction O 2 at z/L = V 1/2 for 
the two situations with Xjj q = 0. 1 and T = 343 K. This figure illustrates the fact that 

the more inert present, the greater is the magnitude of 3X/> /dj. Also, the larger the 

^2 


mole fraction inert, the steeper the oxygen mole fraction gradient is in the z-direction 
at a given J. 


Case D 


As a last example, consider a circumstance where a portion of the water required 
in the current production reaction is evaporated into the gas phase and then recondensed 
in the vicinity of a reaction site. More specifically, for mathematical simplicity it is 
assumed that one-half the total H 2 O stoichiometric requirement evaporates at or near 
z = L, is transported in the negative z-direction, and condenses at a molar rate per 
unit volume equal to that at which the O 2 enters the liquid phase at all z. This is 
equivalent to saying that one-haK the H 2 O reacted at every point arrives in the vicinity 
through the gas phase. 

Thus, qI = iil/4f. 

2 

Assuming the current production is linear with z as in case C, equation (24) yields 

No 

^2 4f6Ae 

and from equation (25) 

NHo = -^^(-f 

^2'-' 4f 6A e V L / 



Substituting these into equation (26) yields for the inert distribution 


X 


I 


Xp exp < 


z 

L 



(33) 
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Equation (27) may now be solved for the reactant distribution 



Both integrals on the right side of equation (34) require numerical solution. 

The first integral on the right side of equation (34) gives the pure gas-phase compo- 
sition. Figure 11 shows Xq for the pure gas phase at 343 K with Xj^ q = 0. 1. The 

2 2 
equivalent result from case C is represented by the dashed lines. 

Consideration of figure 11 yields more information about the situations to which 
cases C and D pertain. Case C requires that the liquid and gas phases be in equilibrium 
with respect to H 2 O at all values of z/L. Thus, the small variation in mole fraction 
H 2 O in the gas phase implies a similarly small variation in liquid-phase concentration. 
If, on the other hand, case D describes the gas phase, the two phases are in equilibrium 
only at z/L = 0, and a very steep concentration gradient exists in the liquid in the region 
of z/L = 1.0. Figure 12 shows the required relation between the mole fraction H 2 O in 
the gas phase and the equilibrium mole fraction associated with the liquid phase (X^ q). 


The second integral in equation (34) gives the additive effect of the presence of an 
impurity. It can be numerically integrated without trouble, but for most purposes it can 
be evaluated directly by making a simplifying approximation involving the exponential 
terms. Combining them, the exponent may be written 


J 



(P 



i i + e 

\ (P 




For J s 0.4 and the 9 - and 0-values listed in table n, the magnitude of this exponent 
is always less than 0. 15. Thus, the exponential may be closely approximated by the 
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quantity 1.0 plus the exponent, and the argument of the second integral thus becomes a 
polynomial in z/L. This is readily integrated analytically and the second term on the 
right side of equation (34) is thus 


6 



k=l 


where 


ai 


= 1 


^4 


= 0 


^2 



a. =J-i(5</^ - 4 - e) 

12 (p 


3-0 — 


1 

3 


slq = - 1 - 0 ) 

® 18 (p 


As in case C, the composition at z/L = ^ 1/2 is chosen as representative. Figure 

13 shows this composition for N 2 as impurity at 343 K. 

The mole fraction O 2 at z/L = f^/2 for pure O 2 and for air is shown in figure 

14 for the same conditions as figure 10. 


Discussion of Cathode Results 

A number of conceivable cathode operating modes have been investigated, this ap- 
proach paralleling that used on the anode. And again, the hope is to have "bracketed" 
most realistically operating cathodes. 

Case A represents a hydrophobic electrode; case B the other extreme with regard 
to wettability and efficient catalyst utilization. Cases C and D embody a linear spatial 
distribution of current, comparing the results with and without H 2 O vapor transport. 

To determine precisely whether H 2 O is transported as vapor - and if so to what 
extent - requires additional analysis of the liquid phase in the cathode. As a first ap- 
proximation, the larger the ratio of cross-sectional area in the gas phase to cross- 
sectional area in the liquid phase, the greater the likelihood for gas-phase transfer of 
H 2 O. Water vapor transport should also be relatively greater at high current densities 
and if considerable amounts of current are produced in the upper (small z/L) regions 
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of the electrode. It is not valid to conclude that vapor transport predominates from the 

fact that Dg q-KOH 0-0 ^ more correct (although still not completely 

2 2 2 

sufficient) comparison is q-KOH 0-0 ’ where C,p is the total 

2 / 2 2 

molar concentration. And C|,D^ q_q /c^D^ q-KOH about 2 at 343 K. 


Again it is found that the parameters J, 6, and cp characterize the gas phase. As 
before, the greater the value of J, the larger the partial pressure gradients in the gas 
phase. Also the smaller the value of cp the greater the propensity of the inert to dis- 
tribute itself equally throughout the gas phase. The parameter 6 is involved only if 
there is vapor-phase transport of HgO, and it may be shown that the greater the value of 
9 the larger the partial pressure gradient of HgO required for a given rate of HgO 
vapor transfer. 

It is worthwhile to note one difference between the gas phase in the anode and cath- 
ode. It can be seen from figures 8 to 10, 13, and 14 that the percent of inerts in the 
cathode will always be higher than that in the gas cavity regardless of the operating mode, 
which is not, in general, true of the anode. This is due to the fact that there is always 
a net flow of gas in the z-direction for low z/L values in the cathode (indeed, for all 
values of z/L unless H 2 O vapor transport is taking place). The effect of this net flow 
is to carry some inerts into the electrode, where they accumulate until the partial pres- 
sure diffusion gradient built up is just sufficient to counteract the net flow effect. 

To characterize, then, any cathode one should first find the value of J. If J is 
less than 0. 1, the mean, or representative, mole percent of O 2 in the electrode will be 
within about 8 mole percent of the gas cavity value, and the assumption of constant re- 
actant partial pressure throughout the gas phase should prove adequate for most purposes 
regardless of inerts or HgO vapor transport. The range 0 s J £ 0. 1 encompasses 
most practical cathodes operating at current densities under 1 ampere per square centi- 
meter. High current densities, closed electrode structures and low operating tempera- 
tures will, however, increase J. 

If an inert is present and J is greater than 0. 1 the values of <p and 9 should be 
calculated. If is small (< 1) a considerable excess of inert in the electrode should 
be expected, and a large value of 9 implies impeded vapor-phase H 2 O transport. If 
these circumstances hold, the mean composition of the gas phase in the electrode can be 
approximated using case C or D or case A if the electrode is known to be highly hydro- 
phobic. 
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SUMMARY OF RESULTS 


A physical model for a fuel- cell electrode has been analysed for various sets of 
operating conditions. By assuming forms for the spatial distribution of interphase mass 
transfer, results have been obtained which yield information of general value concerning 
the composition of the gas phase within the electrode. 

In particular, the three dimensionless parameters J, 9, and cp were found to 
characterize the gas phase in fuel- cell electrodes. The quantity J may be thought of 
as a dimensionless current density, while 6 and q) are ratios of binary diffusion co- 
efficients of the gaseous species present. The complete definition of these parameters 
may be found in appendix A. 

The parameters are useful in estimating the effect of gas-phase conditions on fuel- 
cell performance although they by no means constitute a complete description of fuel- 
cell behavior, since electrical performance depends on a number of factors in addition 
to gas -phase concentration gradients. The maximum or reversible voltage may be cal- 
culated from thermodynamic considerations and depends on operating conditions of tem- 
perature, pressure, and electrolyte concentration. Performance losses, so-called po- 
larization, arise as a consequence of irreversible processes associated with the drawing 
of current from a cell. One of these irreversible processes is the transport of reactants 
and products through the gas phase within the electrodes. Although specific values can 
not be assigned to polarization associated with the gas phase because of the coupled na- 
ture of transport in the cell, it is generally true that the greater the difference between 
the partial pressure of reactant in the bulk gas and that in the electrode, the poorer the 
electrical performance as compared to tl^t theoretically obtainable. And conversely, if 
there is little or no partial pressure difference, performance losses must be attributed 
to other factors such as electrolyte resistivity, kinetic limitations and liquid-phase trans- 
port processes. 

The chief utility, then, of the parameters J, 9 , and cp lies in estimating the 
amount by which the reactant is diluted within the electrode through the setting up of gra- 
dients required to transport reactants and products and by the presence of inert gases. 
And these estimates are useful either as design information or as analytical tools in con- 
siderations of fuel- cell polarization. 

Calculations of J, 6, and cp may be used in conjunction with the following generally 
applicable conclusions which emerge from the various cases considered. 

1. The greater the value of J, the larger the concentration gradients within the elec- 
trode gas phase. Numerical limits have been established for the anode and the cathode 
which allow judgement as to whether gas-phase concentration gradients need be consid- 
ered. More specifically if J ^ 0.05 for the anode, or J ^ 0. 1 for the cathode, the 
maximum variation in reactant mole percent from top to bottom of the gas phase in the 
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electrode will be less than 10 mole percent. 

2. If inert impurities are present, their effect on the gas phase depends on the val- 
ues of 6 and cp, that is, on the nature of the impurity. The larger the value of 0, the 
more transport of H 2 O vapor is impeded. For example, if transport of H 2 O vapor 
takes place, a given rate of water vapor transport requires a larger concentration gra- 
dient for larger magnitudes of 0. 

The smaller the value of cp, the greater the likelihood that the concentration of in- 
ert in the electrode will differ from that in the gas cavity. In the cathode in particular, 
small cp values mean that the percent of inert in the electrode is considerably higher 
than in the gas cavity. 

If quantitative information with regard to inert effects is required, estimates can be 
made from the results presented for the case judged to be most representative of the 
electrode in question. Alternatively, if enough is known about the electrode, the general 
equations can be solved in a manner similar to the cases herein discussed. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 1, 1970, 

120-34. 
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APPENDIX A 


SYMBOLS 
3 


C. concentration of species i, g-mole/cm 


“li 

f 

^vol 




k 

k 

g 

L 

Ni,N. 


n. 


R 


P 

Q 

t 


Xi,x. 


6A 

e 


total concentration, g-mole/cm^ 

2 

diffusivity in binary gas mixture of species i and j, cm /sec 

4 

Farraday’s constant, 9.65x10 A-sec/g-equivalent 
total current produced by electrode volume element, A 

2 

current produced per unit gas-liquid interface area. A/cm 

rly^jL/npjfCrpDj^ jj 6 A (dimensionless) 

2 

2 

permeability, cm /sec 

2 

thermal conductivity of gas mixture, cal/(cm )(sec)(K) 

z- value of greatest penetration of gas phase towards bulk electrolyte, cm 

flux of species i or j in the z-direction, g-mole/(cm )(sec) 

equivalents per gram-mole of reactant 

total pressure, A 

energy flux, cal/(sec)(cm ) 

gas temperature, K 

time, sec 

superficial velocity - volume rate of flow through unit cross-sectional geometric 
area, cm^/sec/cm^ 

mole fraction species i or j (dimensionless) 

coordinate direction normal to electrode plane surfaces; equals 0 at surface 
adjacent to gas cavity, cm 

2 

increment of electrode surface area, cm 

ratio gas-filled cross-sectional area to geometric cross-sectional area (dimen- 
sionless) 


V, 


flux of species i across gas-liquid interface, g-mole/(cm^)(sec) 

Dj_j^/Dj jj Q where subscript R refers to reactant (H 2 or O 2 ) and subscript 
I refers to inert species (dimensionless) 
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P 

T 


viscosity of gas mixture, cP(g/(cm)(sec)) 
interfacial area per unit geometric volume, cm 




(dim ens ionless) 


electrode tortuosity (dimensionless) 
Q_j^ (dimensionless) 

i (1 - 6) — (dimensionless) 

<p 2 


Subscripts: 

I inert 

R reactant 

Superscripts: 
gc gas cavity 

g gas phase 

f liquid phase 
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APPENDIX B 

GAS -PHASE DIFFUSIVITIES 

Binary mass diffusivities were calculated by means of one of two expressions: 
reactant- impurity diffusivities from the Chapman- Enskog kinetic theory formula (ref. 5, 
p. 510); reactant-water- vapor and impurity-water-vapor diffusivities from an expres- 
sion developed by Slattery and Bird (ref. 5, p. 505) from kinetic and corresponding 
states arguments. 

The Chapman- Enskog formula is 


^T^AB 


2,2646X10“^ 

[t 

1/2 


® K “b/J 



2 o 

°"AB^^D,AB 


(Bl) 


where and Mg are the molecular weights of the two species; is one of the 
Leonard- Jones parameters for the mixture and is given by a^g = l/2(a^ AB 

is a dimensionless function of the quantity kT/e^g, where k is the Boltzmann constant 
and is the other Leonard- Jones parameter for the mixture and is given by 


^AB ■ ^^A^B 

(Values of and fig ^g are given in reference 5 (pp. 744-746).) 

The diffusivities predicted by the Chapman- Enskog formula have been found to be 
within 10 percent of measured values for low-pressure (< 10 atm), nonpolar gas mix- 
tures. 

The expression of Slattery and Bird with the ideal-gas law embodied as equation of 
state is 


*^T°AB 


4.43X10"®T 


4/3 

g 


(Pg) (Pg) 

A B 


1/3 


(Tg) (Tp) 

A B 


-3/4 


, 1/2 


M, 


M 


B, 


(B2) 


where (Pp) and (Tp) are the critical pressure and temperature of species i. Val- 

L, i ^ 

ues predicted by this expression are thought to be accurate to within about 10 percent for 
pressures under 10 atmospheres. 
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Numerical values for 0 and cp can be calculated from 


0 = 


^T°I-H20 


<p = 


^T^I-R 


Equations (Bl) and (B2) yield the results given in table III for the species under con- 
sideration in this report. Figure 15 shows as a function of T. In all cases the tem- 
perature is in degrees Kelvin, pressure in atmospheres , and in gram-moles 

per centimeter-second. 
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TABLE I. - VALUES OF 6 
AND <p FOR Hg 
REACTANT^ 



Temperature, K 

i 

343 

473 


e 

<P 

e 

<P 

CH4 

8.0 

0.68 

6. 5 

0.56 

N2 

2.8 

.72 

2.2 

.59 

He 

.68 

I. 5 

.54 

1.21 

Ar 

3. 1 

.78 

2.4 

.61 


and 


TABLE n. - VALUES OF B 
AND q> FOR H2 
REACTANT^ 



Temperature, K 

i 

343 

473 


B 

<P 

e 


CH 4 

2.48' 

0 . 82' 

'2.07' 

0.68 

N2 

.74 

.75 

.61 

.62 

He 

.32 

2.70 

.25 

2. 17 

Ar 

.75 

.72 

.62 

.59 


- ^I-02/®I-H20; 

^ = °I-02/%2°'°2' 
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TABLE m. - EXPRESSIONS TO BE USED DST CALCULATING 6 and ^ 


i 

^T^i-HgO 

^T^i-Hg 

^T°i-02 


195 


Xlg 

0« 

51 0x10" 



2 

CH4 

16,9xlO‘^°T^/^ 

1 5. OxlO' "^T ^/^/n(T/72 . 2) 

5. 33x10' '^T^/^/S7(T/124. 5) 

N2 

51.9X10' 

15. 2xl0“''^T^/^/O(T/59. 0) 
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(b) Sample of electrode interior. 

Figure 1 . - Schematic representation of H2 - O2 fuel cell. 







Mole percent N 2 (Ar) 

Figure 4. - Gas-phase composition near bulk liquid (at Z ■ L) 
case B. Temperature, 343 K; impurity, N 2 (Ar). 
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Figure 8. - Ratio of mole fraotion of inert at 
bottom of volume to .mole fraction of inert at 
top. Case A cathode. 
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Figure 9. - Representative gas-phase composition Ithe com- 
position at z/L=Vi 72). Case C cathode. Impurity, 
N 2 (Ar,CH 4 ); temperature, 343 K. 
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Figure 14. - Mole fraction O 2 at z/L = VT72 
(X^^Q = 0. 1). Case D; temperature, 343 K. 



Figure 15. - Diffusivity function for use in calculating binary gas diffusivities from Chapman-Enskog theory 
(ref, 5, p. 746). 
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